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AND SELF-SIMILAR SOLUTIONS 
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Abstract. We study a nonlinear pseudodifTerential equation describing the 
dynamics of dislocations. The long time asymptotics of solutions is described 
by the self-similar profiles. 



1. Introduction 

In this paper, we study the following initial value problem for the nonlinear and 
nonlocal equation 

(1.1) M( = -|u^|A"m on Rx(0,+oo), 

(1.2) m(x,0) = uo(x) for xeR, 

where the assumptions on the initial datum uo will be precised later. Here, for 
a £ (0,2), A" ~ ((9^/9a;^)"^^ is the pseudodifferential operator defined via the 
Fourier transform 

(1.3) (A^)(o = icr^5(o- 

For a = 1, equation ()1.1|) is an integrated form of a model studied by Head 
[9] for the self-dynamics of a dislocation density represented by Ux- This model is 
a mean field model that has been derived rigorously in [7] as the limit of a system of 
particles in interactions with forces in i. In this model, dislocations can be of two 
types, -|- or — , depending on the sign of their Burgers vector (see the book by Hirth 
and Lothe llj for a physical definition of the Burgers vector). Here, the density 
means the positive density \ux\ of dislocations of type of the sign of Ux- Moreover, 
the occurrence of the absolute value \ux\ in the equation allows the vanishing of 
dislocation particles of opposite sign. In the present paper, we study the general 
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case a e (0, 2) that could be seen as a mean field model of particles with repulsive 
interactions in 

Recall that the operator A" has the Levy-Khintchine integral representation for 
every a £ (0, 2) 

f dz 

(1.4) ~ A"u.(a:) = C{a) / {w{x + z) - w{x) ~ zu;'(x)I{|,|<i} } , 

where C(a) > is a constant. This formula (discussed in, e.g., [6J Th. 1] for func- 
tions w in the Schwartz space) allows us to extend the definition of A" to functions 
which are bounded and sufficiently smooth, however, not necessarily decaying at 
infinity. 

1.1. Main results. First note that equation (|l.ip is invariant under the scaling 

(1.5) u^{x,t) ^uiXx,X°'+h) 

for each A > which means that if u = u{x,t) is a solution to (|l.ip . then u-^ ~ 
u^{x, t) is so. Hence, our first goal is to construct self-similar solutions of equation 
(jl.ip . i.e. solutions which are invariant under the scaling ()1.5p . By a standard 
argument, any self-similar solution should have the following form 

X 



(1-6) Ua(a;,t) = $a(y) with y ^ ^^^^^^^^ , 

where the self-similar profile has to satisfy the following equation 

(1.7) -{a + 1)-^ y^'^{v) = -[K--^a{v))K{v) for all y G R. 

In our first theorem, we construct solutions to equation p.7p . 

Theorem 1.1 (Existence of self-similar profile). Let a G (0,2). There exists 
a nondecreasing function of the regularity C^"*""/^ at each point and analytic 
on {—ya,ya) for some ya > 0, which satisfies 

on (-cx),-yc), 

1 on {ya,+ao), 



(A«$„)(y) = ^ for all ye(-y„,y„). 
a + 1 

Remark 1.2. We can obtain the self-similar solutions corresponding to different 
boundary values at infinity, simply considering for any 7 > and 6 G K. the profiles 
7$Q (^^-i/("+i)y^ + b which are also solutions of equation (|1.7p . 

Remark 1.3. The fact that dy^a has compact support reveals a finite velocity 
propagation of the support of the solution which is typical for solutions the porous 
medium equation, cf. Remark 11.71 below. 
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At least formally, the function is the solution of (|1.7|) . and the self-similar 
function Ua given by (jl.6p is a solution of equation (jl.ip with the initial datum 
being the Heaviside function 



(1.8) uoix) = H{x) 



if a; < 0, 

1 if a: > 0. 



In order to check that Ua given by (|1.6[) solves (|1.1[) , we introduce a suitable notion 
of viscosity solutions to the initial value problem (|l.ip - p.2p . see Section [31 In this 
setting, we show in Theorem 13.71 the existence and the uniqueness of a solution 
for any initial condition mq in BJ7C(M), i.e. the space of bounded and uniformly 
continuous functions on M. Although, the initial datum (|1.8p is not continuous, we 
have the following result. 

Theorem 1.4 (Uniqueness of the self-similar solution). Let a e (0,2). Then the 
function Ua defined in Ij 1.0(1 with the profile constructed in Theorem \l.l\ is the 
unique viscosity solution of equation il.l]} with the initial datum {77 



In Theorem II .41 the uniqueness holds in the sense that if u is another viscosity 
solution to mil), (HHl), then u = u„ on (M x [0, -f oo)) \ {(0, 0)}. 

The self-similar solutions are not only unique, but are also stable in this frame- 
work of viscosity solutions, as the following result shows. 

Theorem 1.5 (Stability of the self-similar solution). Let a E (0, 2). For any initial 
data uq G BUC(K) satisfying 

(1.9) lim uo{x) = and lim ^0(2;) = 1, 

X — ' — 00 X — ^ + 'DC 

let us consider the unique viscosity solution u — u{x,t) of lll.l\) ~ l[77^) and, for each 
A > 1, its resettled version — u^{x,t) given by equation S1.5\) . Then, for any 
compact set K <Z {R x [0, +00)) \ {(0, 0)}, we have 

(1.10) u\x,t)^<^^(^-^j^) in L°°{K) as A ^ +00. 

Theorem 11.51 contains a result on the long time behaviour of solution because, 
first, choosing t = 1 in (ll.lOp and, next, substituting A = t^/^^+i) we obtain the 
convergence of u (a:t^/("+^\t) toward the self-similar profile ^aix). 

On the other hand, convergence (|1.10p can be seen as a stability result when 
we consider initial data which are perturbations of the Heaviside function. This is 
a nonstandard stability result in the framework of discontinuous viscosity solutions. 
It shows that the approach by viscosity solutions is a good one in the sense of 
Hadamard, even if we consider here initial conditions which are perturbations of 
the Heaviside function. 
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Finally, we have the following result of independent interest. 

Theorem 1.6 (Optimal decay estimates). Let a G (0, 1]. For any initial condition 
Uq G BUC{M.) such that uq^x G -^^(I^); the unique viscosity solution u of lll.l]) - nT^) 
satisfies 

||w(-,i)||oo < ||wo||oo and t)||oc < ||uo,x||oo for any t > 0. 

Moreover, for every p G [1, +oo) we have 

(1.11) \M,t)\\p<Cp^Juo^x\\t^ t'^^ for any t > 0, 
with some constant Cp^a > depending only on p and a. 

The decay given in (|l.lip is optimal in the sense that the self-similar solution 

satisfies \\{Ua)x{-,t)\\p = ||(*a)y(-)llp t^^f^. 

Remark 1.7. The equation satisfied by u = li^^ of the following form 

(1.12) vt = {\v\A''~^nv)x 

(with the Hilbert transform denoted by Ti) can be treated as the nonlocal coun- 
terpart of the porous medium equation. Indeed, for a = 2 and for nonnegative v, 
equation (|1.12p reduces to vt = {vvx)x = {^'^/'^)xx ■ ^^^^ ^^'^ porous 

medium equation (see e.g. [23 and the references therein), estimates (|l.lip show 
a regularizing effect created by the equation, even for the anomalous diffusion: if 
vq G L^(M) then v G ^^(M) for each p > 1. Observe also that equation p.l2p has 
the compactly supported self-similar solution v{x,t) = i^"+T<i>^ ^a;/t°+T^, where 
the profile was constructed in Theorem ll.il This function for a — 2 corresponds 
to the well-known Barenblatt-Prattle solution of the porous medium equation. 

Remark 1.8. For a G (1,2), we do not know how to define the product \ux\ (A"u) 
in the sense of distributions, which is an obstacle for us to prove the result of 
Theorem 11.61 in this case, see Section [B] Note, however, that the inequalities from 
Theorem 11.61 are valid for a G (1,2] as well, provided the solution u ~ u{x,t) is 
sufficiently regular. 

1.2. Organization of the paper. In Section [2l we construct explicitly the self- 
similar solution. In Section [31 we recall the necessary material about viscosity 
solutions, which will be used in the remainder of the paper. In Section Hj we 
prove the uniqueness of the self-similar solution. Under the additional assumption 
that the solution is confined between its boundary values at infinity, we prove 
the stability of the self-similar solution, namely Theorem 11.51 In Section [Sj we 
prove further decay properties of a solution with compact support. Applying these 
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estimates, we finish the proof of Tlicorcm ll.Sl in the general case. In Section [6l we 
introduce an e-regularized equation, for which we prove both the global existence 
of a smooth solution and the corresponding gradient estimates. Finally in Section 
[71 we deduce the gradient estimate in the limit case e = 0, namely Theorem ll.6[ 
using the corresponding estimates for the approximate e-problem. 



2. Construction of self-similar solutions 



Proof of Theorem The crucial role in the construction of the self-similar profile 
<i>Q is played by the function 



f2 1) v(x)^l Kia){l-\x\Y' for < 1, 

^ ' ' I for > 1, 

with K{a) = r(l/2) [2"r(l + a/2)r((l + a)/2)]"\ This function (together with 
its multidimensional counterparts) has an important probabilistic interpretation. 
Indeed, if {X{t)}t>o denotes the symmetric a-stable process in M of order a e (0, 2] 
and if r mi{t : |X(t)| > 1} is the first passage time of the process to the exterior 
of the baU {x : \x\ < 1}, Getoor [8] proved that E^(T) = v{x), where denotes 
the expectation under the condition X{0) — x. 

In particular, it was computed in [H Th. 5.2] using a purely analytical argument 
(based on definition p.3p and on properties of the Fourier transform) that A"v £ 
L^R) and 

(2.2) A^vix) = 1 for \x\ < 1. 

Now, for the function v, we define the bounded, nondecreasing, C^+"/^-function 

u{x) = / v{y)dy 



which obviously satisfies u{x) = M{a) for all x > 1 and u(x) — —M{a) for x < ~l 
with 



M{a) = K{a) f\l - \y\T^' dy ^ 
Jo 



2"(a + l)r(i±ii) 
Then, for any (p G C^(R), we can introduce the following duality 

(A"u,^) = / uiy)iA'^ip){y)dy. 
Jr 

This defines A"m as a distribution, because we can check (using the Levy-Khintchine 
formula (|1.4p ) that there exists a constant C > such that 

C\\(p\\w2,'=^m 



|(A»(:r)| < 



1 + \x\^+^ 
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If, moreover, supp ip C (—1, 1), it is easy to check using the properties of the function 
V — v[x) that 

(a,(A"«),^) = -(u,A"(a,^)) = -(u,a,(A"(^)) = (A"(a,u),v?) = (i,^), 

where the last inequaUty is a consequence of (|2.2p . From the symmetry of w, we 
deduce the antisymmetry of u, and then (A"u)(— x) = — (A"u)(.x). Therefore, we 
get the equahty (A"u)(a;) = a; in !?'(— 1,1), however by [12, Cor. 3.1.5], in the 
classical sense for each y e (—1, 1), too. 

Finally, we define the nonnegative function 

1>a(2/) = ^{-(7-^/^"+^^y)+M(«)} with ^-i^MM. 
a+lLV / J a+1 

Now, for Ua — 7i/("+i) = [2^/(0?)]^-'^/*^"+-'^^ we can check easily that is exactly 
as stated in Theorem I LI) which ends the proof. □ 



Let us note that we will not use in the sequel the explicit form of the function 
but only its properties listed in Theorem ll.il 

Remark 2.1. It is known since the work of Head and Louat [TU] (see also that 

1 /2 

the function v{x) ~ K [l — x^) (with a suitably chosen constant K = K{1) > 0) 
is the solution of the equation {A^v){x) = 1 on (—1, 1). This result is a consequence 
of an inversion theorem due to Muskhelishvili, see either [19l p. 251] or [22l Sec. 
4.3]. 

3. Notion of viscosity solutions 

Here, we consider equation (jl.l[) and its vanishing viscosity approximation, i.e. 
the following initial value problem for a e (0, 2) and 77 > 

(3.1) Ut — rjUxx — \ux\ A."u on Mx(0, +00), 

(3.2) u{x,0) = uq{x) for x£R. 

In this section, we present the framework of viscosity solutions to problem (|3.ip - 
(|3.2p . To this end, we recall briefly the necessary material, which can be either 
found in the literature or is essentially a standard adaptation of those results. We 
also refer the reader to Crandall et al. [S] for a classical text on viscosity solutions 
to local (i.e. partial differential) equations. 

Let us first recall the definition of relaxed lower semi- continuous (Isc, for short) 
and upper semi- continuous (use, for short) limits of a family of functions which 
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is locally bounded uniformly with respect to e 

limsup *w'^(a;, i) = limsup u'^(y, s) and liminf i) = liminf u'^{y,s). 

e^O e^O , 

y — >x.s — >t y *t 

If the family consists of a single element, we recognize the use envelope and the Isc 
envelope of a locally bounded function u 

u*{x,t)— limsupu(y, s) and u^{x,t) = liminf u{y,s). 

y — >x^s — yt y ^^1^ *i 

Now, we recall the definition of a viscosity solution for (|3.ip - l|3.2p . Here, the 
difficulty is caused by the measure jzj^-'^^^dz appearing in the Levy-Khintchine 
formula (|1.4p which is singular at the origin and, consequently, the function has to 
be at least C^'^ in space in order that A°'u{-,t) makes sense (especially for a close 
to 2). We refer the reader, for instance, to [20l[3l[T6] for the stationary case, and 
to |151 114] for the evolution equation where this question is discussed in detail. 

Now, we are in a position to define viscosity solutions. 

Definition 3.1 (Viscosity solution/subsolution/supersolution). A bounded use 
(resp. Isc) function it : M x K+ ^ R is a viscosity subsolution (resp. supersolu- 
tion) of equation (|3.ip on M x (0, +oo) if for any point {xo,to) with Iq > 0, any 
r e (0,to), and any test function cj) belonging to C'^(M.x{0, +oo))nL°° {Rx {0, +oo)) 
such that u ~ (j) attains a maximum (resp. minimum) at the point (a;o,io) on the 
cylinder 

Qr{xo, to) R X (^0 - t, to + r), 

we have 

dt(l){xo,to) - ?](l)x^{xQ,to) + \(j)^{xo,to)\ {A"(j){-,to)){xo) <0 (resp. > 0), 

where (A"(/)(-, to))(xo) is given by the Levy-Khintchine formula p.4p . 

We say that u is a viscosity subsolution (resp. super solution) of problem (|3.ip - 
p.2p on R X [0, +oo), if it satisfies moreover at time t = 

w(-,0)<uo (resp. "(•, 0) > (uq)*) • 

A function u : M x R+ ^ R is a viscosity solution of (|3.ip on R x (0, +00) (resp. 
R X [0, +00)) if u* is a viscosity subsolution and is a viscosity supersolution of 
the equation on M x (0, +00) (resp. R x [0, +00)). 

Other equivalent definitions are also natural, see for instance [3]. 

Remark 3.2. Any bounded function u G C^^*^ (with some (3 > max{0, a— 1}) which 
satisfies pointwisely (using the Levy-Khintchine formula (|1.4p ) equation (|3.ip with 
r] = 0, is indeed a viscosity solution. 
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Theorem 3.3 (Comparison principle). Consider a bounded use subsolution u and 
a bounded Ise supersolution v of i3.1\) - [37^) . Ifu{x,0) < uo{x) < v{x,0) for some 
uq e BUC(R), then u <v onRx [0, +oo). 

Proof. Recall that in |141 Th. 5] , the comparison principle is proved for a — I and 
rj — under the additional assumption that uq S W^'°°{M.). Looking at the proof 
of that result, the regularity of the initial data uq is only used to show that 

(3.3) sup ((uo)^(x) — (uo)e(a;)) ^ as e — > 0, 

where {uqY and (ito)e are respectively sup and inf-convolutions. It is easy (and 
classical) to check that p.3p is still true for uq g BUC{M.). The general case can 
be done either considering a variation of the proof of IM] taking into account the 
additional Laplace operator, or applying the "maximum principle" from [16| . or 
following, for instance, the lines of [3j. We skip here the detail of this adaptation. 
This finishes the proof. □ 

Theorem 3.4 (Stability). Let {M^}e>o be a sequence of viscosity subsolutions (resp. 
supers olutions) of equation iS.l]) which are locally bounded, uniformly in e. Then 
u = limsup*u^ (resp. u = liminf, M^j is a subsolution (resp. supersolution) of 
( TO]) on M X (0,+c»). 

Proof. A counterpart of Theorem 13.41 is proved in [3, Th.l]. Here, the result for 
the time dependent problem is again a classical adaptation of that argument, so we 
skip details. □ 

Remark 3.5. One can generalize directly Theorem 13.41 assuming that {u^}e>o are 
solutions to the sequence of equations p.ip with rj = e. Then, in the limit e 0"*", 
we obtain viscosity subsolutions (resp. supersolutions) of equation (jl.ip . We use 
this property in the proof of Theorem 11.61 

Remark 3.6. In Theorem 13. 4| we only claim that the limit It is a supersolution on 
R X (0,+oo), but not on M x [0,+oo). In other words, we do not claim that u 
satisfies the initial condition. Without further properties of the initial data it 
may happen that ?l(-, 0) < Uq is not true. 

Theorem 3.7 (Existence). Consider uq G BUC(R). Then there exists the unique 
bounded continuous viscosity solution u of ^3.1\) - ^i3^) . 

Proof. Applying the argument of |13j (already adapted from the classical argu- 
ments) , we can construct a solution by the Perron method, if we are able to con- 
struct suitable barriers. 
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Case 1: First, assume that mq G M^^'°°(M). Then the following functions 



(3.4) 



u±{x, t) — uq{x) ± Ct 



are barriers for C > large enough (depending on the norm ||uo|jvF2.°°(K))7 and we 
get the existence of solutions by the Perron method. 

Case 2: Let uo G BUOM). For any e > 0, we can regularize mq by a convolution, 
and get a function S VF^'°°(K) which satisfies, moreover. 



Let us call the solution of (j3.ip - (l3.2|) with the initial condition Uq instead of 
uq. Then, from the fact that the equation does not see the constants and from the 
comparison principle (Theorem 13. 3p . we have for any e, (5 > 



Therefore, {u^}e>o is the Cauchy sequence which converges in L°°(]R x [0, +oo)) to 
some continuous function u (because all the functions are continuous). By the 
stability result (Theorem [331), we see that m is a viscosity solution of equation p.ip 
on K X (0, +oo). To recover the initial boundary condition, we simply remark that 
u^{x, 0) — Uq{x) satisfies p.Sp . and then passing to the limit, we get u{x^ 0) — uo{x). 
This shows that m is a viscosity solution of problem p.ip - (|3.2p on M x [0, +oo), and 
ends the proof of Theorem 13.71 □ 

4. Uniqueness and stability of the self-similar solution 

Lemma 4.1 (Comparison with the self-similar solution). Let v he a suhsolution 
( resp. a supersolution ) of equation (OUP with the Heaviside initial datum given in 
m.8\) . Then we have v* < (ua)* (resp. (ua)* < v^,). 

Proof. Using Remark 13.21 and properties of gathered in Theorem 11.11 it is 
straightforward to check that the self-similar solution Ua{x,t) given in (|1.6p is a 
viscosity solution of equation (|3.ip - (|3.2p with the initial condition (jl.Sp . 

Now, we show the inequality (wq)* < w*. Let w be a viscosity supersolution 
of ()3.ip - ()3.2l) with the Heaviside initial datum (II. 8p . Given a > and v°-{x,t) — 
v{a + x,t), we have 



Because of the translation invariancc of the equation (jl.ip . we see that v"' is still 
a supersolution. Moreover, for any a > 0, we can always find an initial condition 
Ua e BUC{R) such that 



(3.5) 



[-"o - wol < £• 



lu" -u'^l <e + S. 



M*{x,0) < (uoTix) < (uoUa + x) < v''{x,0). 



Uq(x,0) < Ua{x) < v"'(x,0). 
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Therefore, applying the comparison principle (Theorem 13. 3|) . we deduce that 

Because this is true for any a > 0, we can take the limit as a — and get (uq), < w*. 

For a subsolution v, we proceed similarly to obtain v* < (wq)*. This finishes the 
proof of Lemma 14.11 □ 



Proof of Theorem \1.4\ We consider a viscosity solution v of equation with the 
Heaviside initial datum (|1.8p . Using the both inequalities of Lemma [4.11 and the 

fact that (wq)* = on (R x [0, +oo))\ {(0, 0)}, we deduce the equality v = Ua 

on (M X [0, +oo)) \ {(0, 0)} , which ends the proof of ThcorcmO □ 

We will now prove the following weaker version of Theorem 11.51 

Theorem 4.2 (Convergence for suitable initial data). The convergence HI. 10]) in 
Theorem \1.5\ holds true under the following additional assumption 

(4.1) lim uo{y) = < uo{x) < 1 lim uo{y). 

y — > — oo y — >+oc 

Proof. Step 1: Limits after resettling of the solution. Consider a solution u of (jl.ip - 
(|1.2|) with an initial condition uq satisfying ()4.1|) . Recall that for any A > 0, the 
rescaled solution is given by u^{x, t) — u{\x, A"+^t). Let us define 

u = limsup*u^ and w = liminf*u^. 

From the stability result (Theorem 13. 4p . we know that u (resp. u) is a subsolution 
(resp. supersolution) of (jl.ip on R x (0, +oo). 

Step 2: The initial condition. We now want to prove that 

(4.2) u{x,0) = u{x,0) = H{x) for .t e R\ {0} , 

where H is the Heaviside function. To this end, we remark that uq satisfies for 
some 7 > the inequality |uo(a;)| < 7 (note that 7=1 under assumption (|4.ip ). 
and for each e > 0, there exists M > such that |uo(a^)| < £ for x < —M. 
In particular, we get 

uq{x) < e + jH{x + M), 
and then from the comparison principle, we deduce 

(4.3) u{x,t) <e + {ul)^{x + M,t) 
with 

(4.4) ul{x,t)^^l[-^^j^) and ^Hy) = [j-'/^"+'^y) . 



NONLINEAR DIFFUSION OF DISLOCATION DENSITY 



11 



Here $2 the self-similar profile solution of (ll.7|l with the boundary conditions 
and 7 at infinity. Moreover, because is continuous off the origin, we can simply 
drop the star * , while we are interested in points different from the origin. This 
implies 

u\x,t)<s + ^l(^ ,i/(.+i) 

and then 

Therefore, for every a; < we have 

Because this is true for every e > 0, we get u{x, 0) < for every a; < 0. We get 
the other inequalities similarly, and finally conclude that (|4.2p is valid. 

Step 3: Initial condition at the origin, using assumption We now make 

use of (|4.f p to identify the initial values of the limits u and u. We deduce from the 
comparison principle that 

0<u{x,Q) <u{x,Q) < 1, 
and then for every a; e M we have 

u{x,0) < H*{x) and u{x,0) > H^x). 
Step 4- Identification of the limits after reseating. From Lemma |4. 11 we obtain 
M < K)* = K)* < u on (R X [0,-|-cx)))\{(0,0)}. 
We have by the construction u <u, hence we infer 

u = u = Ua on (M X [0,+oo))\{(0,0)}. 

Step 5: Conclusion for the convergence. Then for any compact K d {R x 
[0, +oo))\ {(0, 0)}, we can easily deduce that 

sup \u''^{x,t) ~ Ua{x,t)\ ^ as A +00, 
ix,t)eK 

which finishes the proof of Theorem 14.21 □ 

5. Further decay properties and end of the proof of Theorem 11.51 

Theorem 5.1 (Decay of a solution with compact support). Let u be the solution 
to il. l]) - nT^) with the initial datum uq G BUC{M) satisfying for some A> 

(5.1) uo{x) < for |a;| > A. 
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Let also 7 > 6e such that uo(x) < 7 for all a; G M. Then, there exist (3,(3' > 
(depending on a, but independent of A,j) such that 

u{x,t) < Ct-I^, 

and 

u{x, t) <0 for |x| > C't^' 
with some constants C — C{a,A,"f) and C = C'{a, A,"f). 

First, we need the following 

Lemma 5.2 (Decay after the first interaction). Consider o,nd defined in 
Theorem l 1 . li Let v G (1/2, 1) and G (0, ya) he such that ^a{i,v) = v . Let T > 
be defined by 

(^■2) ^l/(a+l)yl/(a + l) ^^''^ 

Then, under the assumptions of Theorem \5.1[ we have 

(5.3) u{x,t)<i^j for all t>T, xeR, 
and 

(5.4) u{x, t) <0 for aU < t < T and \x\> a(i + — 

Proof. Let us denote $2(2/) — l^a (7~^/^"+^^y). Then we have 
^H{x + A)> uq{x) for a; G M, 

where 

^H{x + A)=l\^^^l(^-^±±^J for x + A^Q. 
Now, we apply the comparison principle to deduce that 

*o (^^(^) ^"(^'^) foi" (a:,i)GMx (0,+cx3). 
This argument can be made rigorous, simply, by replacing the function ^H{x + A) 
by (^(2; + ^ + <5)/(iy''"^^'')^ for (5 > and some sequence 0+, and then 

taking the limit 5^0+. 
Therefore we have 

for (x,t)GMx (0,+oo). 
From the properties of the support of we also deduce that 

u{x,t)<Q for x < - (a + y„ {it) ^Z^^+i) ) , 
and then, by symmetry, 

u{x,t)<Q for \x\>A + ya,{-ftf''^"+^\ 
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Moreover, it follows from the monotonicity of that 

7*" ( (^,)il+i) )^"("-^) 
for a; < 0, and by symmetry we can prove the same property for a; > 0. Then for 
T > defined in (|5.2p we easily deduce (|5.3p and (|5.4p . This ends the proof of 
Lemma (lO). □ 



Proof of Theorem \5.1\ We apply recurrently Lemma 15.21 Define = ^7 7o = 7i 
and 

^"+1 = ^» (1 + I;) ' 7n+i='^7n, and -^-—^j^-^=^^. 



This gives 

/4„ = /In ( 1 + 



with 



70 V / I' V 



and therefore 

u" — 1 

w(x,t)<7„ for t >To + ... + r„_i =Xi— -. 
In particular, we get for any n e N 

"(a;, t) < 7oJ^" for t > Kon"- 
with Kq = K/ (/i — 1). This implies 

u{x, t) < 70^0 1"'^ for any t > 0, x£ M, 

with 

Inj^ 

f3^-—>0. 
in /i 

Similarly, we have 

u" — 1 

u{x,t)<0 for |a;|>A„ if t < Tq + ... + T^-i = K- . 

^ — 1 

In particular, we get for any n e N\ {0} 

u{x,t)<0 for \x\ > Ao (^l + , if t<KQi/' 

with = This implies 

u{x,t)<Q for |a;| > ^0(^0)"'^'*^' for ^ > 0, 

with 



In 1 + 1^ ' 

/?' - \ > 0. 

In ^ 

This ends the proof of Theorem 15.11 □ 
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As a corollary, we can now remove assumption ()4.1|1 in Theorem 14 . 21 and complete 
the proof of Theorem 11.51 

Proof of Theorem \1.5[ We simply repeat Step 3 of the proof of Theorem 14. 2i but 
here without assuming (|4.ip . Then, for any e > there exists A> such that 

uq{x) < 1 + £ for |x| > A. 

By Theorem 15 . 1 1 applied to the solution u{x, t) — l — e, this implies that there exists 
a constant C > (depending on e) such that 

u{x,t) < l + e + Ct~f^. 

Therefore, for any for A > the following inequality 

holds true, which implies that u — limsup*u'^ satisfies 

A — *+oo 

u{x,t) <l + e for (a;,t) e M X (0, +oo). 

Since this is true for any e > 0, we deduce that 

u{x,t)<l for (x,t) e M X (0,+oo). 

Let us now define u — min (1, u) . By the construction, 

li{x,t) ^u{x,t) for (a;,i) e M X [0,+oo)\{(0,0)}, 

and, by (|4.2[) . we have u{x, 0) < H*{x) for all x £ R. Therefore, m is a subsolution 
of (|l.ip - (|1.2p on R X [0, +oo) with the initial datum being the Heaviside function. 
Similarly, we can show that u — limsup,u^ satisfies 

A — >-f C30 

u>Q for e R X (0,+oo). 

Hence, the function u — max(0,u), which is a supersolution of (|1.1|) - (|1.2|) on 
R X [0, +oo) with the Heaviside initial datum. 

Finally, the conclusion of the proof is the same as in the proof of Theorem 14.21 
where u (resp. u) is replaced by u (resp. u) . This ends the proof of Theorem ll.51 □ 

6. Approximate equation and gradient estimates 

In this section, in order to prove our gradient estimates of viscosity solutions 
stated in Theorem 11.61 we replace equation (jl.ip by an approximate equation for 
which smooth solutions do exist. Indeed, with e > 0, we consider the following 
initial value problem 

(6.1) Ut = suxx — \ux\^°'u on Mx(0, +oo), 

(6.2) u{x,Q) = uo{x) for a; G R. 
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We have added to this equation an auxihary viscosity term which is stronger than 
A"u and Ux- In the case a E (0, 1], we will see later (in Section[7]) that it is possible 
to pass to the limit e ^ 0+ in L°°(IR), which is the required convergence for the 
framework of viscosity solutions. The difficulty in the case a £ (1,2) comes from 
the fact that, for the limit equation with e = 0, we are not able to give a meaning 
to the product \ux\ (A"w) in the sense of distributions, while it is possible when 
ae (0,1]. 

Our results on qualitative properties of solutions to the regularized problem 
(j6.ip - (|6.2p are stated in the following two theorems. 

Theorem 6.1 (Approximate equation - existence of solutions). Let a G (0, 1] and 

£ > 0. Given (iTiy ifiitidl ddtuui uq G C^(]R.) such that uq^x 

e L^{R)r\L°^{R), there 

exists a unique solution u G C(M x [0,+oo)) n C'^'^{R x (0, +oo)) of i fO]) -ifO |) . 
This solution satisfies 

(6.3) e C([0,T],LP(R))nC((0,T]; VK^'^ (R)) n Ci((0, T], (R)) 
for every p g (1, oo) and each T > 0. 

Theorem 6.2 (Approximate equation - decay estimates). Under the assumptions 
of Theorem \6.1[ the solution u = u{x,t) of I16.1\) - [K^) satisfies 

(6.4) ||u(-,i)||oo < ||wo||oo, \\ux{-,t)\\ 
and 

(6.5) IMMp < CpJ\uo,x\\f^t~'^^''^\ 

for every p £ [l,cx)), all t > 0, and constants Cp,a > (see, \6.2U\) below), inde- 
pendent of e > 0, t> and uq. 

6.1. Existence theory. 

Proof of Theorem \6.1\ Note first that 

(6.6) A"u ^ -A^-^Hux, 

where Ti, denotes the Hilbert transform defined in the Fourier variables by (7iti)(^) = 
i sgn(^) v{£,). We recall that the Hilbert transform is bounded on the i^-space for 
any p £ (1, +oo) (see [211 Ch. 2, Th. 1]), i.e. it satisfies for any function v £ _LP(R) 
the following inequality 

(6.7) \\nv\\,<Cp\\v\\p 
with a constant Cp independent of v. 



16 PIOTR BILER, GRZEGORZ KARCH, AND REGIS MONNEAU 

For a e (0,1), the operator A"~^ defined analogously as in (II. 3p corresponds 
to the convolution with the Ricsz potential A"^^w = Cq,| • * v. Hence, by [lU 
Ch. 5, Th. 1], for any p > 1/a with a E (0, 1] and any function v G L'^{R), we have 

(6.8) \\A°''^v\\p < Cp^a\\v\\q with i = i + l-a. 

Now, if M = u{x,t) is a solution to (|6.ip - (|6.2p . using identity (|6.6p . we write the 
initial value problem for v = Ux 

(6.9) vt = evxx + {\v\A"^^nv)x on Mx(0,+C3o), 

(6.10) v{-, 0) = Wo = uo,x e L\R) n L°°(M) 
as well as its equivalent integral formulation 

(6.11) v{t)^G{st)*vo+ f dxG{s{t ~ t)) * ilvlA^-^Hv) dr, 

Jo 

with the Gauss- Weierstrass kernel G{x,t) — (47rt)^^/^ exp(— a;^/(4t)). 

The next step is completely standard and consists in applying the Banach con- 
traction principle to equation (|6.1ip in a ball in the Banach space 

Xt = C([0,r]; L\R) n L°^\R)) 

endowed with the usual norm ||u||t = supjgjQ 7-](||w(i)|ji + ||w(t)||oo)- Using well 
known estimates of the heat semigroup and inequalities (j6.7p - (|6.8p combined with 
the imbedding L'^(R) n L°°(R) C LP(R) for each p e [l,oo], we obtain a solution 
V — v{x,t) to equation (16. lip in the space Xt provided T > is sufficiently small. 
We refer the reader to, e.g., [II Ej for examples of such a reasoning. 

This solution satisfies (|6.3p for every p e (l,oo) and each T > 0, by standard 
regularity estimates of solutions to parabolic equations. Moreover, following the 
reasoning from IJ, one can show that the solution is regular. 

Finally, this local-in-timc solution can be extended to global-in-time (i.e. for all 
T > 0) because of the estimates ||f(t)||p < Ij^ollp for every p £ [l,oo] being the 
immediate consequence of inequalities (|6.17p . (|6.18p . and (|6.2ip below. □ 

6.2. Gradient estimates. In the proof of the decay estimates of Ux, we shall 
require several properties of the operator A". First, we recall the Nash inequality 
for the operator A". 

Lemma 6.3 (Nash inequality). Let < a. There exists a constant Cn > such 
that 

(6.12) mT^"^ <CN\\A-^MlMr 
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for all functions w satisfying w G L"'^(M) and A"/^w £ L^(R). 

The proof of inequality (|6.12p is given, e.g., in [T71 Lemma 2.2]. 
Om' next tool is the, so called, Stroock-Varopoulos inequality. 

Lemma 6.4 (Stroock-Varopoulos inequality). Let < a < 2. For every p > 1, we 
have 

(6.13) / iA"w)\w\P~^wdx> 1^^^ f U'^lwl^y dx 

for all w £ LP{R) such that A'^w e LP{R). // A"w e L^{R), we obtain 



(6.14) / (A"u;) sgnw dx > 0. 
Moreover, ifw, K"w £ L^(]R), it follows that 

(6.15) ( (A"w)u;+ dx > and ( (A"u;)w" dx > 0, 
Jr Js. 

where — max{0, w} and — max{0, —w}. 

Inequality (|6.13p is well known in the theory of sub-Markovian operators and its 
statement and the proof is given, e.g., in J18', Th. 2.1 combined with the Beurling- 
Deny condition (1.7)]. Inequality ()6.14|) . called the (generalized) Kato inequality, is 
used, e.g., in |t6j to construct entropy solutions of conservation laws with a Levy dif- 
fusion. It can be easily deduced from [6l Lemma 1] by an approximation argument. 
The proof of (|6.15p can be found, for example, in [TSl Prop. 1.6]. 

Remark 6.5. Remark that inequality (|6.14p appears to be a limit case of (|6.13p 
for p — 1. Inequality (|6.15p for follows easily from (|6.14p . by a comparison 
argument, if for instance w G Cjj"(R). Finally, remark that the constant appearing 
in (|6.13p is the same as for the Laplace operator d^/dx^ — — A^. 

Our proof of the decay of v{t) — Ux{t) is based on the following Gagliardo- 
Nirenberg type inequality 

Lemma 6.6 (Gagliardo-Nirenberg type inequality). Assume that p £ (l,oo) and 
a>Q are fixed and arbitrary. For all v G L^{R) such that A"/2|^;|(p+i)/2 ^ L'^{R), 
the following inequality is valid 



(6.16) \\v\\l < Cn A"/2|z;|(p+i)/2 

where 



V 



2 



1; 



p(p + a) , pa + 1 

a = , — , 

p — 1 P ^ 1 

and Cn is the constant from the Nash inequality 116. 12\) . 
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Proof. Without loss of generality, we can assume that ||v||i ^ 0. Substituting 
w — |i,|(p+i)/2 in the Nash inequality (|6.12[) we obtain 



|(j9+1)(1 + q) 



2 



"ll(p+l)/2- 



Next, it suffices to apply two particular cases of the Holder inequality 

< \\v\U, as weU as ||«||(,+i)/2 < Mf'^+'^My^'+'K 



V 



1 



and compute carefully all the exponents which appear on the both sides of the 
resulting inequality. □ 



Proof of Theorem \6.SX The first inequality in (16. 4|) is the immediate consequence of 
the comparison principle from Theorem l3.31 because classical solutions are viscosity 
solutions, as well. Maximum principle and an argument based on inequalities (|6.15p 
(cf. [13 for more detail) lead to the second inequality in (|6.4p . We also discuss this 
inequality in Remark 16.71 below . 
For the proof of the L^-cstimate 

(6.17) llu.Wlli < ||mo,.||i 

(i.e. (|6.5p with p = 1 and Cp,a = 1), we multiply equation (j6.9p by sgnu = sgnu^: 
and we integrate with respect to x to obtain 

^ [ \v\dx^e ( v.^^sgnvdx+ [ {{A°'''^nv)\v\)^sgnvdx. 

The first term on the right hand side is nonpositive by the Kato inequality (i.e. 
(|6.14p with a = 2) hence we skip it. Remark that (formally) 

{{A"-^nv)\v\)^sgnvdx ^ [ {A"-^nv)v^{sgnv f + {A''-^Hv^)v dx 



{{A°'-'^nv)v) ^ dx = 0. 

Now, approximating the sign function in a standard way by sgns{z) — zj^J + 5, 
integrating by parts, and passing to the limit (5^0"'", one can show rigorously that 
the second term on right hand side of the above inequality is nonpositive. This 
completes the proof of (|6.17p with p = 1. 

Next, we multiply equation in (|6.9p by |ti|P^'^i; with p > 1 to get 
l_d 
V dt 

We drop the first term on the right hand side, because it is nonpositive by ()6.13p 
with a = 2. Integrating by parts and using the elementary identity 



\v\^dx^e \ v^:,|f|P^2„dx+ / {{A''-^Uv)\v'^^\vY-'^vdx 
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we transform the second quantity on the right hand side as follows 

{{A''-^Hv)\v\) \v\P-^vdx = -^^ f iA"v)\v\P-\dx. 

Consequently, by the Stroock-Varopoulos inequality (j6.13p (with the exponent p 
replaced hy p + 1), we obtain 

(6.18) lll«WII^^-fT^|K"''l^l^^"^^'t- 

Hence, the interpolation inequality (|6.16|) combined with ()6.17p lead to the following 
inequality for ||w(i)||p 

(6.19) ^^Hm; < (c.ii.„ii^+^)/(-^))"^ i\\vm;Y''-"'^^'-'' ■ 

Recall now that if a nonnegative (sufficiently smooth function) f — f{t) satisfies, 
for aU t > 0, the inequality f'{t) < -Kf{t)^ with constants A' > and /3 > 1, 
then f{t) < {K{P - l)^)"^/^''"^). Applying this simple result to the differential 
inequality (I6.19p . we complete the proof of the L^-decay estimate ()6.5p with the 
constant 



(6.20) = C 



,_i4p(a+l) 



where Cat is the constant from the Nash inequality (|6.12p . □ 

Remark 6.7. Note that, for every fixed a, we have limp_>oo Cp_Q = +oo. By this 
reason, we are not allowed to pass directly to the limit p — s- +oo in inequalities 
(j6.5p (as was done in, e.g., [TTl Th. 2.3]) in order to obtain a decay estimate of v{t) 
in the L°°-norm. Nevertheless, using (|6.19p we immediately deduce the inequality 
||''^('; t)\\p ^ ll'^o(')llp valid for every p e (1, oo). Hence, passing to the limit p +oo 
we get 

(6.21) \H-,t)\\oo<\\vo{-)\\oo. 

In general, we cannot hope to get a decay estimate of ||w(., t)|joo better than that 
in (|6.2ip . because each constant is a solution of equation (jl.ip . Moreover, one can 
show that if v is constant on an interval at the initial time, then it will stay equal to 
the same constant on an interval depending on t, because of the finite propagation 
phenomenon that can be seen for the self-similar profile. 

7. Passage to the limit and proof of Theorem 11.61 

Now, we are in a position to complete the proof of the gradient estimates (jl.lip . 
First, we show that form the sequence {u^}e>Q of solutions to the approximate 
problem (|6.ip - (|6.2p one can extract, via the Ascoli-Arzela theorem, a subsequence 
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converging uniformly. Theorem [33] on the stabihty and Remark 13.51 imply that the 
limit function is a viscosity solution to (|l.ip - (|1.2p . Passing to the limit e ^ 0+ in 
inequalities (|6.4p and (|6.5p we complete our reasoning. 

Proof of Theorem[LE First, let us suppose that uq e C°°(K) n W'^'°^{R) with 
""0,2; G i^(R) n i°°(K). Denote by u'^ = u'^{x,t) the corresponding solution to the 
approximate problem with e > 0. 

Step 1: Modulus of continuity in space. Under this additional assumption, we 
have 

(7.1) ||<(.,t)||p<Cpt-> 

with Cp = Cp,q||wo,2;||i'°^^' and 7p = ( 1 — M. The Sobolev imbedding theo- 
rem implies that there exist some (3 G (0, 1) and Co > such that 

(7.2) \u'{x + h,t)-u'{x,t)\ < IhfCoCpt-'^''. 



Step 2: Modulus of continuity in time. Let us consider a nonnegative function 
E C°°(M) with supp(p C [—1, 1] such that Jj^ (p{x) dx = 1, and for any S > set 

ips{x) — S^^ip{S^^x). Then, multiplying (|6.ip by ips and integrating in space, we 

get 

_d 

dt 



u^(a;, t)ips{x) dx ] = e(u^(-, t), {(ps)xx) 

ips{x) \ul{x,t)\{H A""i ul{x,t))dx, 



and then with l/p+ l/p' = 1 
d 

(7.3) 



— ( / u^{x,t)ips{x)dx 



<s\\u%;t)\U\{ips)xx\\l 

+ ||^,-||oo||<(-,t)||p||||H A"-l Ul{;t)\\p,. 



Here, we have used relation (|6.6p . Combining inequalities (|6.7p and (|6.8p with 
estimate (|7.ip , we get for p' > 1/a 

\\n A"-i <(-,t)||p, < Cp,Cp,,c.Cgt-^-. 

Then for any bounded time interval / C (0, +c») there exists a constant Cj^s such 
that for all t E I, we have for any e £ (0, 1] 



< \s\Ci,s- 



Now, for any t,t + s E I, we get 

u^{x,t + s)lps{x) dx — / u'^{x,t)ips{x) dx 
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Therefore, the following estimate 

\u'{0,t + s)-u''{0,t)\ 

<\s\Ci^s+ / Vs{x)dx- sup {\u%x,t + s) -u%0,t + s)\ + \u^x,t) -u'{0,t)\) 

holds true. Using the Holder estimate (|7.2p . we deduce that there exists a constant 
C/ depending on /, but independent of S and of e e (0, 1], such that 

\u'{0, t + s)- u'{0, t)\ < \s\Ci,s + CiSf^. 

Since the above inequality is true for any 6, this shows the existence of a modulus 
of continuity luj satisfying 

\u''{0,t + s) -u^O,t)\ <uJi{\s\) for any t,t + sel. 

By the translation invariance of the problem, this estimate is indeed true for any 
X e R, i.e. 

(7.4) \u%x,t + s) -u%x,t)\ <uji{\s\) for any t,t + sel, xeR. 

Step 3: Convergence as e — > 0+. From estimates (|7.2p and (|7.4p . and using 
the Ascoli-Arzela theorem and the Cantor diagonal argument, we deduce that 
there exists a subsequence (still denoted {u^jg) which converges to a limit u g 
C(R X (0, +00)). By the stability result in Theorem 13.41 (see also Remark l3.5p . we 
have that u is a viscosity solution of (|l.ip on K x (0, +00). 

Step 4-' Checking the initial conditions for uq smooth. Remark that for uq e 
^2,00 ^^j^ i-j^g barriers given in p.4p with some constant C > uniform in 
e G (0, 1]. This ensures that u is continuous up to < = and satisfies u{-, 0) — uq, 
so this proves the result under additional assumptions. 

Step 5: General case. The proof in the case of less regular initial conditions 
simply follows by an approximation argument as was in the proof of Theorem 13. 71 

Step 6: Cradient estimates. To pass to the limit e — > 0+ in estimates (|6.5p . we 
use the inequality 

(7.5) h-'\\u'{- + h,t) - u'{;t)\\p < \\uU;t)\\p 

with fixed h > 0. Hence, by the Fatou lemma combined with the pointwise conver- 
gence of toward u, we deduce from (|7.5p and (|6.5p that 

h-^\\u{- + h,t)-u(-,t)\\p<Cp,o,\M\f^t-^('"-p) 

for all h > 0. For every fixed t > 0, the sequence {h^^{u(- + h,t) — u(-,i))}/i>o 
is bounded in LP{M.) and converges (up to a subsequence) weakly in LP{R) toward 
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Ux{-,t) (see, e.g., |2li Ch. V, Prop. 3]). Using the well-known property of a weak 
convergence in Banach spaces we conclude 



hx(-,i)||p < liminf + h,t)- u{-,t)\\p < Cp,c«ho,,||f^t-^(^^i). 
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